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ABSTRACT 


The  problem  of  the  approximation  of  a  given  random  signal  with 
modulated  pure  tone  signals  is  studied  from  the  standpoint  of  approximation 
of  the  povrer  spectrum.  The  equal  energy  and  moment  matching  techniques  are 
developed  for  the  selections  of  the  modulation  parameters  of  the  pure  tone 
sources.  Techniques  are  developed  for  the  selection  of  parameters  of  a 
svfitched  pure  tone  source.  The  power  spectinim  is  also  approximated  by  a 
rational  polynomial  scheme,  which  uses  a  mlnimax  procedure,  developed  in  the 
appendices,  for  selection  of  pure  tone  operating  parameters.  Several 
examples  of  each  of  these  techniques  are  given. 
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DTERODOCTIOI 


IQ  a  prerioua  report  [9]  the  theoretical  prohleas  aaaoclated 
vlth  the  approxlJMtlng  of  random  acoustical  signals  vlth  pure 
tone  sirens  has  heen  treated.  This  report  vlth  the  more 

practical  aspects  of  the  engineering  appiroxlaatlon  of  random 
signals.  A  nuBd>er  of  teehnlq]aes  tore  deTeloped  and  applied  to 
a  particular  signal  to  Illustrate  the  methods. 


Manuscript  released  "bj  author,  Octoher  1964,  for  publication  as 
an  BTD  Techni cal  Documentary  Bepart. 
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I.  POWER  SFBCTULL  iraSITT  IBBHICnBlTIOW 


GlTen  a  contlnuoua  record  of  finite  length,  it  is  not 
possible  to  estlnste  the  sutocorarlanee  function  C(t)  for 
arbitrary  t.  Thus  in  place  of 


C(t) 


11a 


1 

*4 


T/2 


z(t)x(t  -¥  T)dt 


(1.1) 


ve  nay  calculate 


1 

^oo^^^  *  T  -|t|  J  )  3c(t  +  5’)dt  (1*2) 

n 


vhere  is  known  as  the  apparent  aatoeorarlance  function. 

Here  T  is  the  length  of  the  record  which  is  restricted  by 

|t|<i^<t  ,t  being  the  aaxiauB  Isg  tine  to  be  consld- 

ered.  VSeal  ^th  the  problm  of  the  relation  of  T  to 

T  in  section  IH.  * 

n 

Let  ])(e)  be  a  function  defined  by 

(o  |t|  >  T 
m 

low  define 

C^(t)  .  d(t)  ejt) 

It  is  clear  that  C  (t)  is  calculable  for  a  glren  rsoerd, 
and  edso  that  C  Is  defined  for  |f  |  >  T  ,  even  though 
^(t)  was  not  defined  in  this  region. 

If  we  assume  that  the  random  process  is  ergodlc  then 
the  average  of  a  large  nueber  of  finite  records  is  equiva¬ 
lent  to  a  single  average  taken  over  a  single  record  of 
infinite  length.  Thus  ve  can  write: 

ave(C^(T))  -  D(t)C(t) 

It  follows  then  that  this  relation  will  have  a  well  defined 
Fourier  Integral  of  the  form 
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ETC  (P^j(f))  -  Q(f)  *  P(f)  P(i-f)di 

vhere  Q(f)  !•  tb»  Foutlor  transfozn  of  D(t)  and  P(f)  is 
ths  Fourier  transfoza  of  C(t)  ,  texned  tbs  power  spectral 
densl'ty.  The  sysibol  denotes  tbs  convolution.^  The 
Fburler  transfom  0(f)  of  a  function  g(i)  is  defined,  as 
usual,  h]r 

0(f)  .*^di 

irtiere  f  ■  2irw  .  Thus, 

ave  {P^(f))  -^Q(f-l)P(l)di.  (1.3) 

We  then  identify  the  average  of  P  (f)  as  a  smoothing 
(average  over  frequency)  of  the  ^ue  power  spectral  densl-ty 
with  a  weight  function  Q(f-i). 

We  now  turn  our  attention  to  the  problem  of  non>contlnuous 
records  of  finite  length.  Suppose  that  we  have  a  function  x(t) 
recorded  on  a  record  of  finite  length  with  the  function  z(t) 
specified  only  at  a  number  of  equally  spaced  points  on  the 
time  axis  t  .  Let  us  also  asstmie  that  x(t)  is  given  only  over 
equally  spaced  Intervals  of  difference  At,  l.e. 

t  *  0,At,  2At, . .  .nAt 
then  C(t)  can  only  be  estimated  for 
|t|  *0,  At,2At,...nAt 


If  E(f)  is  the  Fourier  transform  of  h(t)  and  a(f)  is  the 
Fourier  transform  of  g(t),  then  the  Fourier  transform  of 
Mt)g(t)  Is  given  as 

0(i)E(i>f)di  -  y^0(i-f)E(l)di 
(1)  P*  for  a  further  discussion. 


G(f)*E(f)  -jT 
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Tb«  integral  eqimtlon 


2P^(f)  coe  (2irft)df 


|t|  ■  q  ■  0,1, ...n. 


il.k) 


If  soluble  at  all,  mist  be  satisfied  by  a  function  P.(f) 
vhleh  vanishes  for  f  >  f  ■  ,  even  thou^pi  the  original 

power  spectral  density  extends  beyond  |f^|  .  "Wii*  Intro¬ 
duces  a  new  problma,  that  was  not  encountered  In  the  case 
of  continuous  spectra,  naaely  that  of  aliasing.^  Consider 
the  two  slgMLls  shown  In  Figure  1. 


Fig.  1.  -  SasqpUng  of  two  slnusodal  waves 

We  note  that  In  tews  of  the  diacrete  saeple  x(t),  t  >  .2q, 
q  >  0,1,2,...,  the  two  waves  of  different  frequency  both  give 
the  sew  result.  Thus  we  see  that  equally  spaced  tine  ssaples  of 
one  sine  wave  could  have  cobm  frost  any  of  aaay  other  sine  waves. 
Special  care  aust  therefore  be  exercised  in  regard  to  statistical 
properties  of  a  continuous  function  x(t).  In  particular,  the 
relation  of  the  actual  power  spectral  density  P(f)  and  the  one 


^  [2],  p.  35 


k 


i^ch  WB  are  capable  of  laeaaurlng,  '  xepraaent  a  critical 

problem. 


To  pxerent  such  a  dlfflcul'ty,  ve  vlU  aaaume  that  our  samp¬ 
ling  time  is  so  chosen  that  the  maximum  freq]aene3r  component  of 
x(t)  is  less  than  \/2£A  i.e.,  ' 

1^1 

This  then  forces: 

r  Pa(') 


It  should  be  remarked  that  this  is  seldom  the  ease  in  actual 
practice. 

One  would  now  ask  -  "Is  there  a  better  sampling  scheme  than 
that  of  equally  spaced  samples?"  If  one  Indeed  exists^  it  appears 
that  it  has  escaped  attention.  When  the  problem  of 
becomes  important,  it  is  usual  to  follow  some  method  of  filtering 
and/or  smoothing  the  data  to  moke  appropriate  modifications  in 
the  computations.  Ve  make  no  attempt  to  do  this  in  the  present 
wozk.  References  [2],  [3]  and  [^]  are  of  Interest  in  approaching 
this  problem. 


|f|  <  (2At)‘^ 
|f|  >  (2At)"^ 


n.  IBWERICAL  IglKBWIMkTIOB  OF  POWER  SraCTBA 


Let  us  denote  the  walues  of  the  sample  function  detexnlned 

at  times  t  "^t,  q  ■  0,1,2, . .  .n  as  x  ,  that  is  we  define 

% 

-  x(q^t)  ,  q  -  0,1, . .  .n 

It  is  easily  seen  that  the  expression  analogous  to  (1.2)  is  now 


2  C(iAt) 


n-r 


X  X 


q  qBr 


q.-0 


where  tut  >  x^t  and  rM>,l,2, . .  .a  <  n  . 


(2.1) 


To  dateneine  the  power  spectrum  then  we  take  the  appropriate 

analogue  to  the  relation  (l.V)  and  find  that 

n-1 


r  •  0,1,2,...^ 


The  fhctor  of  2  appeara  in  e«i.  (2.2)  dne  to  the  ftet  ttet  «« 
vleh  to  cooelder  only  thoee  power  epeetra  STanetrle  about  w  -  0. 
Thue,  glren  a  particular  record,  ve  hare  a  maeerlcal  aethod  to 
ccapute  the  corarlaace  function  and  the  poear  apectral  denalty. 

If  there  are  aereral  dlatlnct  recorda  to  he  uaed,  aa  la 
uaually  the  came,  «e  define 


C  -  are  (C  ) 

*  It 


F 

r 


are  (P  1 

T 


(•.3) 


where  the  areraoe  la  taken  for  each  r  orer  each  dlatlnct  record. 

Fy  the  property  of  ergodlclty,  which  la  tacitly  aaawMd,  we  are 

aaaxned  that  la  the  Uiat  of  Ineroaalnely  larse  nvebera  of 

recorda  that  will  alnoat  aarliMy  conrerge  to  C(t)  at  polnta 

of  definition  of  7  ,  For  alapUclty  ffon  here  on,  we  drop  the 

hara  on  C  and  f 
r  r 


It  la  alao  q^te  frequently  dealrahle  to  know  the  prohahlU'^ 
denal^  for  a  glren  randon  function.  The  probability  deiMl^ 
function  F(y)  la  defined  aa 


P(y)<iy  -  Proh(y  <  x  <  y  ♦  dy|. 


Ita  detemlnatlon  for  our  diacrete  aai^Ung  then  tnma  out  to 
he  nerely  a  couatltB  prooeaa. 


m.  PIAMMUtt  FOR  TBB  MglSCRBglB 


Blackean  and  TUkey  [2]  hare  approached  the  prohlen  of  the 
aaount  of  data  required  about  a  particular  prooeaa  In  order  to 
■ake  atatoMnta  concerning  the  reltabniy  of  eatlaatea  of  power 
apectra  and  corarlance  functlona.  Becanae  of  the  aatrononleal 
nwber  of  calculatlona  larolred  In  eatlmtlag  a  power  apectnM 
®  dealred  preclalon.  It  la  laportant  that  we  coaqprteilae 
between  reaolu^^on  and  atablU^.  We  ahall  only  reproduce 
certain  parte  of  their  reaulta.  The  reader  will  find  the 

oF  arrlTlng  at  theae  reaulta  In  auction  A.  23  A.  28, 

B.  23 -»B.  2%  of  [2]. 
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W»  define  ee  laeed  In  eectlon  I,  •• 


(reeolutlon  In  e.p.e.)  «  ~  «r 

T 


(3.1) 


lAere  ie  Meeured  In  eeeoade,  "Iteeolatlcm’'  le  a  aeaeuxe 
of  the  eoneentxatlon  of  a  epeetral  eetlaate  In  frequency  unite, 
l.e.  a  aeaeure  of  the  vldth  of  a  frequency  In  ubleh  no 
atteapt  le  aade  to  further  dleerlnlnate.  We  eaa  e:q>reee  the 
etahlUty  aeeoelated  vlth  a  particular  eetlnate  In  texM  of  a 
epread  In  db,  db  ■  20  log^.  Interral  containing, 

nth  preecrlbed  probrt>lll^,  ^  »tio  of  true  power  to  eeti- 
aated  power.  To  thle  end  we  defend 


^  ^  2.*  230 _ 

(dOjt  Bunge  In  d^ 


k  -  1  -f 


boo 


(  90^  Bunge  In  db) 

k 

*  (9^  Buigo  In  db) 


(3.2) 


k  -  1  4> 


840 


(989^  Bunge  In  db) 


^y  90  per  cent  range,  it  met  be  e^phueized  that  we  hare 
9  out  of  10  chaneea  of  finding  each  indlTldual  eetlnate 
nthln  the  preecrlbed  db  range  of  Ite  average  value. 

The  nwdier  of  recorde,  P,  of  length  T  required  to  yield 
the  particular  atabillty  and  reeolutlm  choaen  le 
ezpreeelble  ae: 

I  (1  +  k) 

n  ■  3 

We  find  that  we  need 


(3.3) 


5T  f 
n  • 


lage  for  the  coBqputatlon  of  the  covariance  and  At  -  l/Jf 

aa  the  else  of  the  time  Intervale,  f  la  the  i^Trtimim 
frequency  component  preaent  in  x(t)  .“*fcing  thie,  we  can 
conpute  the  required  nonber  of  data  polnte  to  be: 

n  -  (1.5  k  +  P)T  f 


(3.4) 


We  now  give  eeveral  illuetratlve  ezaaqplee. 
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1. 

SuppoM  that  v»  art  gXrmn  raeorda  of  langth  1  aaeond  Iwrinc 
fraq^ancy  eonpoaenta  19  to  10^  e.p.a.  aad  daalza  a  raaolutlon  of 
30  e.p.a.  We  vlah  to  knov  bov  aany  aaaplea  of  thla  racord  anat 
he  Bade  to  hara  a  eonfidenee  leral  of  90  per  cent  at  +  2  db. 

Thua  aa  ealeolate: 

k  -  101 

P  -  1 

B  ■  600 

n  -  3  X  10^ 
va  then  need  only  1  racord. 

BxMople  2  • 

S^9Poae  aa  are  glTan  recorda  of  len^^th  2  aecoada.  We 
wlah  to  find  out  boa  nany  plecea  Ta  alll  need  to  bare  a  98 
per  cent  confidence  at  ^  1  db,  altb  frequency 

eonponenta  of  10  e.p.a.  and  altb  50  e.p.a.  reaolntlon. 

T  •  2  X  lO"® 

B 

k  -  8V1 

P  -  1^.2  ^  1» 

4'  recorda  are  then  raq;nlred. 

17.  AMPLITULB  AID  raBftUJaUi  MOIOIAIBD  Stabata 

We  hare  at  our  dlapoaal  K(  <  25)  loa  fratoency  and  L(<  9) 
hl^  freq^iency  aln^  tone  ganeratora  (alrena)  altb  ahlcb  aa 
vlah  to  ^n?i^3clinate  a  given  nunerlcal  poaar  apectral  denal'^. 

airan  nay  be  either  aapUtude  or  fre^piency  nodulated 
althin  certain  reatrlctlva  raoKea. 

let  ua  find  the  poaar  apectral  denaltlea  for  each  of  the 
three  nodea  of  operation  of  a  alngle  airan: 
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a)  Pure  tone:  X^(t)  coa 

(1^.1) 

then  p^(«)  «  &(»-\) 

Is  the  center  frequency  of  the  siren,  and  A^  Its  pover  and 
&(<e-«i^)  Is  the  usijal  Dirac  delta  function  defined  as 


such  that 


jT  ft(w-«^)d»  ■  1 


h)  Aoplitude  Modulation 


then 


:^(t)  COS  O^t)  cos  u^t 


\M  -  -j—  ♦  a^)  +  6(»-^-a^)) 

♦  \  6(— 

•fET  i»  called  tlM  modulation  motor  and  (W2lr  the 
modulation  frequency. 

c)  Frequency  ItidBiatlcm 

Vt)  cos  («^t  +  sin  Q^t)  ; 


then 


+  Jq  (\)  8(«-«j^)) 


(^.3) 


^  is  tezmed  the  deviation  ratio,  and  ft /2rr  the  modulation 
frequency.  J  (6)  Is  the  Bessel  functlim  of  order  n  and 
argument  5. 

Figures  2  and  3  give  representative  pictures  of  power  spsotra 
of  the  pure  tone  and  amplitude  modulated  signals. 
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P(U)) 


I  i 


A 


k 


(jJ 


Fig.  2  -  Power  apeetrum  of  a  pure  tone  signal 


Fig.  3  -  Power  spectrua  of  an  aoplitude  nodulated  signal 

The  power  spectrum  of  the  freq>iency  modulated  spectrum  Is  not 
as  slnqply  Illustrated.  In  Figure  4  we  Illustrate  seweral  oases 
of  frequency  modulation.  A  table  glrlng  the  absolute  magnitudes 
of  each  frequency  ccmponent  Is  presented  In  appendix  C. 

Sxippose  we  are  glren  a  particular  power  spectral  density, 

3,  that  is  specified  numerically  at  each  of  the  points  w  •  ndM. 
Let 

S  S(n/:ea)  ;  n  ■  0,1,2,... 
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<kir  obJactlTs  iM  to  approzlnate  S  }fj  tho  use  of  E  lov 
fxoqiiioiicy  slrBiio,  l.«.  opoxmting  In  a  ftrafoency  range  of  0 
to  |/2t  and  L  high  frequency  alrena,  operating  in  a  range 
|/2ir  to  I  ,  each  operating  In  one  of  the  }  a»dee  glren  ehore 
by  a),  b)f  or  e).  Each  of  the  three  typea  of  Modulation 
abore  depend  upon  the  epeelfleation  of  a  center  frequency 
w  .  To  detenlne  w  ve  sake  the  requirenent  that  each  of  the 
K  ,  and  einilarly  the  L,  elrene  hare  the  mum  total  power 
output.  Ihla  aaeuaption  is  wade  In  order  to  assure  that 
any  criteria  used  for  the  selection  of  paraMters  will  lead 
to  realisable  power  distributions  aaong  the  sereral  sirens. 
The  total  power  to  be  distributed  aaong  the  K  low  frequency 
sirens  is  given  by 


P 


L 


/■ 


a 

Sta)dM 


S 

4 


V  m 


and  slallarly  for  the  L  high  frequency  sirens 

T  ^ 


vs  bare  assuoed  that  8(m)  is  essentially  sero  for  w  >  {  . 
Suppose  that  Q  and  are  respectively  the  operating 
powers  of  the  low  and  nigh  frequency  sirens.  Then 


u. 


I.l.lll. 

8  -  5.5 


I  .  .  I 

6  -  .1 


..lll.l.llni 

5  - 


5-1.  5-6. 


5-8. 


.  il  I  I  I  .  1 1 .  I  ■  li  I  . 1 1  .1 1  I  I  .  . 

5  -  10. 


5-2.5  5  -  15. 


8-3. 


I  I  .  I  .  1 .  1 .  1 1 


5-20. 


...  1 1 

L 

1 1 1 1  i  1 » 1 1  ■  1  ■  1 

ii 

1... 

5-25. 

Fig.  k  -  Pover  spectrum  of  a  frequency  modulated  signal  for  various  5 
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I?  Qg  is  aloore  tbs  aUovmble  physical  output  of  a  high 
frsqusncy  siren,  then  it  is  necessary  to  scale  S  in  such 
a  wanner  that  is  reduced  to  this  limiting  saliie.  Ve 
similarly  scale  8  if  is  outside  the  allosshle  operating 
range  of  the  lev  frequency  sirens.  The  Interval  (0,$  )  nov 
needs  to  be  broken  up  into  K  equal  energy  intervals*  i#t 
<V  \+  **  *^  f««»nC3r  lnt.r«l  inrtich  th. 

sirens  operates.  Then  is  determined  by: 


subject  to  the  conditions  that 


0 


\+l  "  ^ 

Similarly  se  section  the  high  frequency  range  by  requiring 


vhere 


1  *  K^l, ...,  I« 


^K+1  “  ^ 
^K+IH-1  “  ^ 


Thus  ve  have  partitioned  the  power  spectnsi  into  parts 
(see  Figure  ^). 
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Fig.  ^  •  Energy  pertltlonlng 


Let  us  nov  examine  one  frequency  tend,  es  lUuetreted 
in  Figure  6. 


Fig.  6  -  A  single  frequency  bend 


4 


W6  nov  tlM  res'trlc'tloii  'tbft't  only  ono  slxvn  Is  sUowed 
to  opersts  In  each  band. 

We  quote  a  powerful  theorem  of  fourler  analysis; 


on  ^e  Interval  [a,b],  then  the  sequence  of  functions  • 

)  C  T  tends  uniformly  to  X(w)  as  n  -» «•  at  each  point  w  c 
n  n  — 


of  continuity  of  Xe  C  ^  the  Fourier  coefficient  defined  by 
— *  21  ^ 


C  -  /  T  (t)X(t)dt 
n  J  n 


The  ISchebyscheff  polynomials  foxm  such  a  system  of  functions. 
Then  n*^  Tchebyscheff  polynomial  is  a  polynomial  of  degree  n 
In  the  free  variable,  that  is 


n 


Th-ttig  this  fact  we  then  can  claim,  by  the  properties  of  unlfozm 
convergence,  that  If  we  force  all  the  moments  of  a  function 
3  to  be  the  same  as  all  the  moments  of  another  function  3  then 
3^  and  3  are  the  same  functions.  Guided  by  this  we  Intend  to 
\ise  the  moment  matching  between  the  actual  power  spectrum  and 
the  assumed  power  spectrum  as  the  criteria  for  selecting  all 
the  parameters  that  need  to  be  deteralned.  Generally  we  vlll 
not  have  enough  free  paqmmeters  In  order  to  force  all  of  the 
moments  to  be  equivalent.  In  this  case,  we  will  arbitrarily 
choose  the  lowest  monents  needed  to  specify  the  parameters. 

let  us  apply  this  scheme  to  the  amplitude  modulated  signal 
given  by  e^.  (4.2).  The  pure  tone  signal  vlll  then  be  the 
sane  except  A  is  set  identically  zero.  Writing  the  balance 
of  the  zeroth^  first,  second,  and  fouirth  moments  gives  us: 


gae  any  book  on  advanced  mathematics  such  as  [1].  It  turns 
out  that  we  must  restrict  x(t)  such  that  t\(t)  is  square 
Integrable  for  each  m.  However,  we  need  not  worry  about  thla, 
since  our  x(t)  vlll  be  bounded  and  continuous  on  a  finite 
Interval. 
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•^i+l 


P,/K  1  <  K 
K 


S  s(w)dw  ■  A  +  ■  Q  ^  “ 

f,  i  2  i  jp^L  1> 

/  MB(u)dii>  ■  (a  +  ^  )  u 

J  1  2  1 


r^i+1 

J  ) 


s(»)d«  . 


/' 


1+1  _  k 

S(«)d»  - 


We  can  then  write  the  parametere  A  ,  A  ,  a  ,  «  aa 

1'  l'  i'  i  * 


A 

w  ■  — 

1  Q. 


0. 


*  X 

A 

Ira 


(>*.5) 


A.  ■ 


iHiere: 


5  -  B  -  2i»^A  + 

T  ■  D  -  W.C  +  6«  ^  -  W  ^A  + 

*  w  *  1  i  1 

A  - 

B  ■  (A«)^y*2 
*1 
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•2 

C  m  (£4#)^y 

**1 


D  ■ 


H- 


Sl 


Sl±l 


■til 

Th\is  we  hawe  detenilned  the  parameter*  of  the  1  alren  given 
that  we  vleh  to  approzinate  the  i>ower  •pectrum  with  an  amplitude 
modulated  algnal. 


In  the  case  of  frequency  modulation,  the  situation  Is  not  quite 
as  stral^tfhrward*  vs  again  have  four  parameters  to  select, 
namely,  A^,  w^,  A^,  5^  .  First  notice  that  the  only  place  that  5^ 
occurs  Is  In  the  argument  of  the  Bessel  function  J  *  It  Is  fairly 
ohvlous  that  there  will  he  no  slq>le  method  of  choosing  6  hy  the 
solution  of  the  moment  matching  equations*  we  must  choose 

somewhat  arbitrarily.  First  notice  that  equation  for 

does  not  change  when  referred  to  the  frequency  rather 
amplitude  modulation  cases.  Thus  Is  readily  calculated.  By 
making  reference  to  Figure  k  we  txy^to  choose  the  5  idilch  will 
most  closely  give  the  same  kind  of  gross  characteristics  as  the 
actual  power  spectral  densl-ty  exhibits  over  the  (f  ,  f  )  Interval, 
low  that  5^ has  been  selected  we  proceed  In  a  similar  to  the 

method  used  In  the  derivation  of  relations  (4.^).  Thus  we  get 


Aw 


i^-7) 


A  -  — Si— 

idiere  we  have  used  the  notation  of  equations  (^.6)  and 
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^]r  xise  of  the  relation 


2  2  2  ^ 

*  ■  ♦  *2  ■  ^1^2  ® 


ve  can  ahov  that 


J°(8)  -  1 


^(8)  -i 


nraa  throuj^  equations  (M>*7)  vs  can  determine  the  paraaeten 
associated  vlth  a  frequency  nodulated  signal  to  qn>TOxlnate 
given  power  spectrun  on  an  Interval 

As  can  be  seen  by  cos^arlng  Figure  4  with  Figure  the  use 
of  frequency  modulation  Is  capable  of  giving  nucn  wider  distributions 
of  the  power  over  toe  Interval. 

▼.  BAirooMij  swucHKD  mao. 

Consider  the  form  of  the  power  spectrum  of  a  pure  tone  siren, 

l.e. 


If  this  signal  is  turned  on  and  off  In  a  random  fashion  we  will 
have 


cos  w,  t 
k 


(5.1) 


where  (^(t)  Is  a  random  function  assvming  values  4-1  or  *1,  as 
Illustrated  In  Figure  7* 


Fig.  7  ->  Bandon  on/off  ampUtude  modulation  function 

We  denote  the  rwan  aifitchlng  time>  l.e_^>  the  expected  time 
between  changes  in  value,  of  Qj^(t)  by  If  we  assume  that 

the  period  between  switchings,  i.e.  t  in  Figure  7>  is  a  Poisson 
dlstrubuted  random  variable,  then  ^  a  strai^tfprwMCd  calcu¬ 
lation,  we  find  that  the  power  spectral  density  of  the  signal 
(!}.l)  is  given  as 


where 


3^(-) 


) 


Hote  that  this  is  the  first  time  that  we  have  been  lead  to  a 
non-dlsciete  power  spectral  density*  This  non-discrete  form 
is  realised  due  to  the  random  nature  of  Q*  (t) .  In  previous 
methods  we  have  always  been  concerned  with  deterministic  forms 
(see  section  IV)  idilch  necessarily  lead  to  discrete  spectra. 

Since  we  are  only  concerned  with  power  spectra  that  are  syumetrlc 
about  w  «  0  ,  we  write  P(w)  as 


P(«)  -  I  Sj^(«) 


2ir 


<1 .  ) 


2 


(5.5) 
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Blis  amounts  to  a  folding  about  u  >  0  of  the  spectrum.  Figure  8 
shoes  a  typical  form  for  P(w)  as  given  by  (5 •5) 


Fig.  8  -  Power  si>ectrum  of  a  randomly  svltched  pure  tftiv*  signal 

Bie  half  height  width  of  the  curve,  as  shown  In  Figure  8,  Is 
given  as 

Oils  shows  that  the  width  of  the  spectrum  Increases  and  the 
height  decreases  as  the  mean  switching  •M"**  decreases . 

One  could  now  set  up  a  least  square  method  for  the 
approximation  of  a  given  si>ectrum  by  a  system  of  on-off 
modulated  sirens.  This,  however,  would  be  a  rather  involired 
scheme,  since  each  siren  Involves  3  parameters  («  ,  T  ,  A.  ) 
and  we  are  lead  to  a  system  of  nonlinear  equations  for  their 
determination.  When  we  examine  a  particular  spectrum,  there 
will  quite  probably  exist  certain  resonances  that  are  similar 
In  appearance  to  the  one  Illustrated  In  Fig.  9.  Suppose  Fig. 

9a  represents  the  spectrum  we  are  Interested  In  approximating. 

We  wish  to  remove  the  resonance  by  approximating  this  portion 
of  the  power  spectrum  by  an  on-off  modulated  signal.  We  go 
through  the  following  steps  to  find  the  values  of  the  parametezer 

(l)  Draw  a  straight  line.  A,  approximating  the  background 
under  the  peak. 


(2)  Determine  the  maximum,  6  , 
value  of  P  at  the  point  1. 
the  peak  on  the  line  A. 


value  obtained  and  6^,  the 
B  is  roughly  hal  f  way  wider 
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P<co) 


(b)  -  Modified  apectrum  ahovlng  the  removal  of  the  reaonance 
by  meana  of  an  on/off  modulated  alren 

Fig.  9 

(3)  Drav  a  line,  D,  through  the  point  C,  parallel  to  the 

line  A.  The  point  C  la  determined  by  P  ■  )/2 

and  placed  half  vay  between  the  aidea  of  the^pett. 

(4)  Meaaure  the  width,  r,  of  the  peak  on  the  C,  and 
redetemlne  the  point  C. to  be  in  the  center  of  the  peak. 

(3)  Betexmlne  the  parametera  of  the  alren  aa 


a  w  of  the  point  C  (3*^) 

2Tr 
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By  MSigalng  th«  parameters  as  given  In  Eqp.  to  one  of 

the  sizans  ve  are  left  vlth  a  spectrum  as  Illustrated  in 
Figure  This  reduced  spectrum  is  then  more  likely  to  be 

easily  approximated  by  some  other  method. 

The  same  procedure  vith  a  little  more  freedom  of  choice 
can  also  be  used  for  two  overlapping  peaks  idiexe  the  maxima 
are  distinct;  but  close  together. 

VI.  QEHERAL  RKPRESHTmiQB  ^  m  BATOOM  SICMAL  m  APPBOgUMATICMI 
^  m  PPHER  SPECTBOM  n[  RATIOBAL  PQUHCICIALB 

In  a  report  by  Samuels  [9]  the  problem  of  representing  a 
random  signal  by  a  series  of  orthogonal  functions  vlth  ortho¬ 
gonal  random  variables  as  coefficients  has  been  studied.  Here 
ve  shall  be  concerned  vlth  the  solution  presented  there  vhen 
the  power  spectral  densi'^,  S(w),  has  been  approximated  by 


a  ratio  of  polynomials  as  foUoM: 


The  first  computational  problem  associated  vlth  finding 
a  solution  is  the  approximation  of  the  actual  power  spectral 
density  by  a  form  similar  to  eq,\iation  (6.1).  Ve  have  two 
major  alternative  methods  by  idilch  ve  can  attack  this  problma: 

1)  We  can  use  a  least  square  procedure  in  vhich  ve  determine 
the  constant  ^  minimizing  the  Integral 


vlth  respect  to  each  of  the  a^  and  b^  .  v(w)  is  the  velf^tlng 
function  and  3(w)  is  the  actual  power  spectral  density. 
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2)  We  can  uae  a  alnlnax  procedure  in  vhlch  the  coefficients 
a^  and  are  determined  by  minimising  the  difference  betseen 

the  actual  function  and  its  approximation. 


The  least  squares  method  is  presented  in  detail  in  Appendix  1, 
vhile  the  Mlnimax  procedure  is  presented  in  Appendix  2.  While  the 
least  squares  method  gives  an  over  all  measure  of  the  fit  of  the 
approximating  functlcm,  one  is  not  assured  that  very  large  errors 
vlU  not  be  present  in  localised  areas.  The  Mlnimax  procedure 
remedies  the  objectionable  character  of  the  least  sq]uares  method, 
in  that  ve  are  assured  that  it  vlll  not  miss  local  structure  by 
large  amounts.  However,  ve  loose  the  possibill'ty  of  approximating 
the  function  in  a  ^best**  sense  over  the  total  region  with  the 
particular  solution  generated.  Thus  a  decision  must  be  as 
to  idilch  characteristic  must  be  more  faithfully  reproduced: 

a)  local  characteristics  idiile  possibly  lasing  overall  fitting; 
or  b)  gross  characteristics  while  possibly  losing  local  effects. 


Once  the  form  of  equation  (6.1)  has  been  deteimlned,  l.e.  a^ 
and  b^  are  known,  we  can  detcimlne  the  orthogonal  functions  for 
the  approximation  of  the  random  signal.  The  cooqputational  scheme 
that  is  used  to  determine  the  parameter  of  this  representation  is 
given  by  the  following  six  steps: 

a)  Solve  the  polynomial 


for  the  roots  Q  .  Mote  that  there  vlll  be  exactly  p  roots  of 

vhlch  p/2  of  these  vlll  be  the  negatives  of  the  other  p/2  roots. 

Denote  these  independent  roots  as  a  ,  k  >  1,2,... p/2  .  It  must 

also  be  recognized  that  ve  can  not  Iror  the  present  determine  (l 

except  as  a  function  of  X  . 

n 

b)  Determine  the  roots  of 


These  roots  viU  come  in  pairs  (one  the  negative  of  the  other). 
Let  us  denote  the  roots  for  vhlch  the  imaginary  part  is  positive 
by  k  ■  1,...  p/2  .  The  values  of  vlll  not  be  dependent 
on  X^,  thus  they  are  readily  competed. 


c)  Find  the  quantities 
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k  *  1^2^ • • •p/s 


d)  Solve  the  nonlinear  system 


a  -1» 

nl  1 


0  ,-i«o 

nl  2 


“ni-n/a 


0  --i«T 

nl  1 


a  -+i«, 

nl  1 


fl  ,+iWT 

nl  1 


-1 


-1 


a  ,+i«- 

ni  1 


a  AiM 
ni  1 


-V‘“l 


(6.3) 

«  0 


for  X  ,  vbere  ve  have  substituted  the  expressicms  obtained  for 
a  infection  a).  T^ie  is  the  equation  50  of  [9]»  Number  the 
solutions  according  to  the  Index  n.  Since  (6.3)  vill  be  a  non> 
Hn«»A-r  equation,  there  may  possibly  be  an  infinity  of  solutions. 


e)  Solve  the  systm  of  linear  equations. 
p/2  + 


l,...p/2 


j>+  (flnk+i«8)T 

)  ®sfe - 

BjmX  nk  n 


^  Aaki 


-  (-anJi+i^)T 


-a  . +i« 

nk  n 


.)  -  0 


S  ■  !,...» p/2 


for  each  A  given  each  of  the  X  determined  in  section 

d).  Note  ^2feat  one**of  the  unJaunnas^is  not  detemined.  let  this 

one  be  ^  ^nk'^^x^/2'  ^  ^ 

and  A^A^^yg,  k-l,...p/2  is  known  for  each  X^. 
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f)  Oliren  all  the  foregoing  ee  then  have 
p/2 

♦  (t)  .  y  Ul* 

n  Z_i  nk  nk  ^ 


kvl 


By  forcing  to  be  imitary,  l.e.. 


I 


f(t)  ♦  (t)dt  -  1 
n  n 


ve  detexmlne  A 


‘np/2  . 

Ihus  <^.4)  representB  the  systaa  of  orthonormal  functions 
necessary  to  approximate  the  given  random  signal  x(t)  as 


-Z 

nso 

irtiere  6^  are  orthogonal  random  variables. 


It  takes  little  examination  to  see  that  this  will  require 
a  very  lengthy  computation  for  reasonably  sized  p  and  q. 

To  Illustrate  the  application  of  the  method  of  rational 
polynomial  approximation  in  obtaining  the  orthogonal  function, 
we  consider  the  simplest  case  possible,  namely:  p  e  1  and  q  -  0  . 
Ihus  S(w)  is  approximated  by 


correlation  function 
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This  Is  tbs  power  spectral  densi^  of  a  process  with  a  Maricoff 
autocorrelation  function  Illustrated  In  Fig.  10.  After  some 
manipulation  (see  [9]  p.  17-20)  we  find  that 


w  cos(^^)  +  p  sin  ■  0 

n  2  2 


Thus  the  X  are  detexmlned  and  ue  can  vrlte  t  as 
n  ^n 


"  Allir  »  t) 

n  n  n  n 

n 

irtiere  a  Is  determined  as 
n 


cos  w 
n 


t  4*  ^  sin  w  tl 
n  • 


(6.5) 


Tn.  APPUCATIOa  TO  A  PARTICOIAll  SIOIAL 

As  an  illustration  of  the  methods  presented  in  the  previous 
sections,  we  now  give  several  applications  to  a  particular  signal. 

Figure  11  illustrates  the  representative  trace  of  the  random 
signal  that  has  been  digitalized  and  used  in  this  analysis.  From 
this  signal  we  have  calculated  its  amplitude  prohahillty  function, 
see  Fig.  12,  its  covariance,  see  Fig.  1$  and  its  power  spectral 
density,  see  Fig.  14,  according  to  the  methods  of  section  II. 

It  is  to  be  noted  that  the  power  spectral  density  ie  quite 
ill -behaved.  This  is  due  to  the  nature  of  the  random  noise 
generator  used.  In  practice  one  would  expect  a  less  harsh 
cxirve.  For  some  of  the  analysis  we  have  used  a  "Smoothed" 
power  spectrum  which  would  probably  be  closer  to  the  types  of 
spectra  physically  measured.  This  is  exhibited  in  Figure  1^. 

Note  that  there  exist  marked  resonances,  however,  they  are  not 
as  pronoxmced.  Also  the  rapid  low  level  backgrovmd  has  been 
filtered  out.  When  this  power  spectrum  is  used  it  will  be 
clearly  noted. 
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A.  yrequency  and  Aaplltode  Modulation 


Aa  in  our  first  illustration,  let  us  suppose  that  ve  vish 
to  approxiinate  the  powar  spectral  densi^  vith  the  pure  tone 
sirens  only.  13ien  according  to  the  equal  energy  partitioning 
and  the  moment  matching  technique  of  section  17 ue  can  easily 
compute  the  operating  frequencies  and  powers  of  each  siren. 
Again  we  wish  to  point  out  that  the  equal  energy  method  of 
assigning  these  parameters  to  each  unit  leads  to  realizable 
power  assignments  to  each  unit  and  in  addition  offers  an 
Intuitively  rational  method  of  selection.  It  is  rational  in 
that  one  of  the  foraoost  considerations  of  approximating  a 
random  signal  is  that  the  energies  associated  with  particular 
frequency  regions  be  as  faithfully  reproduced  as  possible.  This 
assures  that  the  response  characteristics  of  any  system  which 
uses  the  approximate  signal  as  an  input  is  as  close  to  its 
response  when  s\d>Jeeted  to  the  physically  occurring  signal. 

One  of  the  questions  raised  is:  how  many  pure  tone  sirens 
are  necessary  to  simulate  the  signal.  As  an  illustration,  a 
number  of  eases  have  been  carried  out  on  the  power  spectrum  of 
Figure  14.  The  following  list  gives  the  number  of  pure  tone 
sources,  the  table  on  which  the  parameters  am  exhibited,  and 
the  appropriate  figure  cm  idilch  a  comparison  is  made  between 
the  approximation,  the  actual  or  smoothed  power  spectra. 


» 


Case  HO. 

Ho. 

of  sirens 

Table 

Figt 

Hi^ 

Low 

1 

25 

9 

1 

16 

2 

20 

6 

2 

17 

3 

15 

6 

5 

18 

4 

10 

6 

4 

19 

5 

5 

4 

5 

20 

*See  Appendix  D. 

Along  the  MSseissi.  of  the  Figures  l6  throu^  20  the  vertical 
lines  represent  the  center  frequencies  at  which  the  sirens 
function.  The  straic^t  line  segments  represent  the  power  of 
the  siren  smoothed  over  its  Interval  of  operation,  and  the 
irregular  line  indicates  the  real  power  spectrum.  It  is  clear 
from  Figures  l£  and  17,  that  cases  1  and  2  approximate  the 
actual  spectrum  quite  closely.  Cases  3  and  4  altho\::gh  not 
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•pproxlnatliig  the  actual  pouar  apeetrua,  do  rapreaant  the 
■Booth  power  •pectrua  rather  veil.  The  last  case  falls  to 
give  anything  other  than  an  Indication  of  the  gross  power 
dlstrlhntlon  and  Is  cosq^letely  devoid  of  any  local  structure. 

As  another  aeasure  of  the  aeanlngfuUness  of  these  approzlaatlonSf 
the  aapUtude  prohahlUty  densl'ty  has  teen  computed  for  ease  1, 
see  Figure  22.  When  cosq^ared  to  that  of  the  original  signal, 
we  find  that  the  fit  Is  q,nlte  good  except  for  the  very  hl(^ 
aaplltudes.  The  tables  associated  vlth  ease  1  throu^  ^  also 
give  the  appropriate  frequency  and  amplitude  modulation  parameters 
for  the  Intervals.  Ve  can  use  these  parameters  to  distribute  the 
energy  over  the  Interval  In  a  manner  that  vlU  more  closely  yield 
the  approximate  structure  given  by  the  figures.  It  should  be  clear 
through  these  examples  that  the  number  of  sirens  necessary  Is 
directly  a  function  of  the  power  spectrum  that  no  general  rule 
can  be  stated. 

B.  Bandomly  Switched  Pure  Tone  3^gn^^ 

For  the  present  consider  the  smoothed  power  spectrum  of 
Figure  1^.  Note  that  there  are  five  marked  resonances  Indicated 
by  the  arrows.  Each  of  these  has  structure  similar  to  that  of 
Figure  9a.  We  wish  to  remove  these  resonances  by  use  of  randomly 
■witched  pure  tone  sirens. 

After  following  the  procedure  outlined  In  section  T,  ve 
arrive  at  the  following  parameters  for  the  five  sirens  required 
to  remove  these  resonances: 


Resonance 


No.  k 

\ 

’k 

1 

3 

289.0 

1.000 

2 

IB 

314.1 

.500 

5 

26.5 

351.8 

.571 

4 

48 

326.7 

.500 

5 

59.5 

395.8 

.222 

When  the  parts  of  the  spectrum  are  resioved  which  are  contributed 
by  these  sirens,  we  have  left  the  spectrum  represented  In  Figure  22 
by  the  lower  line,  termed  the  resolved  spectrum.  The  resolved 
spectrum.  Figure  23,  Is  much  smoother  will  be  far  easier  to 
approximate  with  other  methods  than  the  spectrum  ve  started  with. 
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C.  Rational  Polynomial  Fitting  «nd  Tlite  of  tb»  Maricoff  Correlation* 

Rjr  removing  the  zeaonanoes  In  the  aaootheat  power  spectrum  hy 
the  methods  of  the  prevloas  subsection^  we  have  arrived  at  the 
resolved  spectrum^  Figure  23*  Oils  has  the  same  shape  as 

that  of  Figure  10  and  thus  suggests  that  It  he  approximated  hy 
the  use  of  an  orthogonal  expansion  In  the  sense  of  section  VI. 

We  find  that  a  curve  given  hy 


P(«) 


48600 

6o  <f  w 


(7.1) 


fdlls  Just  helov  the  resolved  spectrum  (see  Figure  23) .  v  choose 
such  a  curve  so  that  the  portion  left  after  subtraction  vUl  he 
everywhere  positive*  Oie  left-over  can  th«»n  he  approximated  hy 
other  means.  Ohlng  the  length  of  record  analysed,  10  sec,  ve 
find  that  eigenfrequency  equation  (6*5)  Is  in  the  form 


2Tr 

60  "n 


+  tan  (lOiru  )  ■  0 
n 


(7.2) 


vhere  Is  now  in  cycles  per  second.  It  is  clear  that  the 

roots  of  this  equation  are  hounded  as 


<  IOtt  u  <  ^iH-2) 
2  n  2 

or 


—  <  ^  ^  — 

20  n  20  10 

•ttus  the  frequencies  will  he  spaced  at  20  per  cycle.  points 

out  the  haslc  problem  with  the  approximation  of  the  actual  signal 
hy  the  orthogonal  function  schaie  of  section  71.  gamely, 
that  it  leads  to  xunrleldy  nvad>ers  of  sirens  operating  within  the 
frequency  regions  of  Interest* 

This  same  power  spectrum,  namely  that  of  equation  (7.1)  can 

he  approximated  hy  the  use  of  only  one  siren  ^dilch  Is  switched 

on  and  off  in  the  maimer  of  section  V,  if  the  parameters  are 
chosen  as 
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w  0 
T  50 


SIOOtt 


Remembering  that  the  signal  of  the  randomly  switched  siren  is 
of  the  form 

x(t)  -  —  [IfQ(t)]  cos  wt, 
then  the  approximating  siren  has  the  modulation 

i(t) 

For  physical  operation^  u  may  haTe  to  be  chosen  to  be  small, 
but  not  zero.  If  this  amplitude  is  outside  the  physical 
operating  limits  of  one  siren,  then  several  sirens  may  be 
used  to  generate  the  amplitude,  idiile  all  are  similarly 
modulated  with  (1  +  Q(t)).  This  type  of  operation  also  has 
its  drawbacks  in  that  it  is  doidttful  that  the  very  moitn 

switching  time  can  be  physically  realized. 
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Fig.  11  laical  time  track  of  the  functloQ  uaad  In  Section  VH 


.10 


Fig.  12  AinpUtude  probal>llity  density  for  tbs  Bampled 
signal 
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C(at) 


13  Autocorrelation  function  for  the  senqpled  signal 


280r 


Fig.  14  Power  spectrml  deneity  of  the  sampled  signal 


Fig.  Sinoothed  power  spectxml  densl'^ 


P(U)) 
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Fig.  l£  Pover  spectral  density  approxlaatlon  vlth  34  pure- 


280  r 
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17  Povtr  •peetrml  density  epproxLMitloin  vltb  26  puz««tone  •< 


SMOOTHED  POWER  SPECTRUM 


o 

o 

CVJ 


o 

lO 


o 

(O 


& 

o 


o 

o 


38 


Fig.  l3  Pomr  spectrAl  density  approzlmatlon  vlth  21 
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Fig.  19  Power  epectrml  density  approodLoation  with  l6  pure-tone  sources 
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Power  spoctral  deiielty  approximation  vith  9  pure^tone 
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Appendix  A 


laaet  Sy^'****  Deteimlnatlon  of  Rntlonfcl  Polynf  1>Te  Coefflclente 

T|ie  problem  vhlch  ve  vleh  to  consider  le  that  of  approximating 
a  given  tabulated  function  f(x)  In  tezna  of  a  ratio  of  two  poly¬ 
nomials  In  X,  i.e.,  ve  construct  a- function, 


such  that  f  (x)  Is  the  ’'best”  (in  some  specified  sense)  approximation 
of  f(x).  " 

In  this  appendix  ve  vlsh  to  examine  the  problem  of  the  determin¬ 
ation  of  the  coefficients  under  the  so-called  "least  squares” 

criterion.  As  visually  formulated,  the  least  squares  method  calls 
for  the  minimisation  of  the  function 


Q(a,p) 


tf(t)  -  yt)l^ 


¥(t)dt 


vlth  respect  to  the  parameters  fit,  g,  ,  idierc  V(t)  Is  soew  non-negative 
veightlng  fanctlon  idilch  many  authors  set  arbitrarily  to  unity.  Form- 
using  our  svqpposed  f  (t)  ve  have 


■/' 


[f(x) 


l^w(x)dx 


(A.l) 


^or  ^  ^  minimised  vlth  respect  to  the  parameters  a^, 

it  Is  necessary  that 


1  ai  0,1,. ..p 


^1 


-  0 


1  *  0,1, . . .q 

kk 


(ik,2) 


Of  course  the  eolution  for  (A.2)  nuet  he  farther  esamlned  to 
Mcertaln  vhlch  solution  sctually  giTes  the  nlnlinm.  Perfom- 
Ing  the  indicated  differentiation  on  Q/o^)  as  giren  in  (A.l) 
it  is  seen  that  it  is  desirable  to  choose  a  selecting  function 
V(x)  such  that  se  can  remo^  the  unkxiosns  ^  fron  the  denoninator. 
Pjr  this  considerAtion  se  are  1*4  to  oske  tie  choice  of  V(x)  as 

"(i)  - 1  y  (A.j) 


i^re  K  is  an  aihitraxy  positire  predetexained  constant.  ¥e 
propose  to  call  this  choice  the  na-tn-y*  sainted  least 

squares.  The  reason  for  including  the  tem  e'^^wlU  be 
apparent  later.  Clearly,  W(x)  as  given  in  (A.J)  is  non- 
negative.  Q(g^gj  can  then  be  written  as: 


QCflbfiJ  -  J  [f(x)^p^x^ 


-K5t^ 
e  dx 


¥e  first  observe  that  if  f(x)  is  bounded  as  x  -*  •  ,  then  we 
must  require  that  p  <  q  for  ^.4)  to  be  bounded.  Also  se 
note  that  the  system  of  paraiasters  (g^  can  be  detexnlned 
only  within  an  arbitrary  multiple.  Obus  without  loss  la 
generality  we  set 


(A.4) 


a  >  1 
o 


We  have 


or 


P 


'Z  V^>  r(x)xV“ax 


(A.5) 


and 


e  dx 


or 


.  i+J  -lit- 
x)x  e  4x 


■L< 


I 


xV*^ax 


(A  .6) 


for 


J  -  1,2,... p 


Vov  define 


'nfrl 


<  ^  f^(x)  x%  *^dx 

*D.l 


and  recall  that 


(A.7) 


I 


n  -KX 
X  e  dx  • 


Dfl 


-Ky 

It  is  nov  clear  shy  se  chose  the  tem  e  in  the  vei^tlng 
function  (A. 3)*  Using  the  notation  of  (A.7)>  (A.6)  and  (A.3) 
can  be  rewritten  in  the  form: 


P 


L  'L 


i+j+i  "  ^j+i 


J  ■  0,1, . . .p 


J  -  1,2,... P 


(A.8) 


k6 


Qjr  defining 


K* 

, 

i+J*P 


1  ■  •  •  pp  }  J**  •  •  pp 

1  •  P^l> • • pP^I^ 1  i  J  *  ♦ • pp 

1  B  1, ...p  ;  J  B  pfl, ...pfqfl 

1  B  •  •  •P^^^X  P  j  B 


j  B  X^  0  0  •P 
J  B  pfX,  .*.pfqt^ 

J  B  X/  •  .  .P 

J  B  pfX, ...pfq+X 


ve  can  then  vrite  (A. 8)  in  matrix  notation  as  the  set  of 
Xlnear  equations: 


pfqtX 

I  V.-. 


i  •  e  •  e 


(A.9) 


Note  that  if  ve  had  not  chosen  the  velf^tlng  function  as 
ve  did;  then  ve  vouXd  not  have  heen  l«d  to  a  Uno^'r  system  of 
equations  for  the  unknowns  $  ,  but  to  a  nftTiHw«»B'r  system.  For 
ntimerlcaX  work  ve  prefer  a  Xlnear  system  to  a  nonXlnear  one. 


Appendix  B 

PeTelopment  of  the  Mlnlmax  Procedure 

Before  development  of  the  actual  procedure  for  minimitv 
approximation  of  arbitrary  functions  vlth  rational  polyncmlals, 
ve  first  give  several  theorems  and  results  to  place  the  procedure 
on  a  firm  mathematical  basis. ^ 

let  It  be  required  to  approximate  a  given  function  g,  vhere 
the  real  values  y^  •  g(t^)  at  the  H  distinct  points  tj^  for 
1*0,  ...  H-l  constitute  the  only  Infoxmatlon  given  conceming 
8* 


let  the  set  T  consist  of  the  H  points  t^^,  1  *0,  ...  N-1, 
irtiere  tj^  <  t^^j^  .  !Bie  Index  set  S  with  M  elements  vlU  be 
used  to  indicate  a  subset  of  T  consisting  of  the  distinct 
points  t,  irtiere  1-0,  ...  M-1  and 


2^1 


<  <  »-l 


thus 


‘o  ^  ‘•1  ■=  ".1.1  ^  Vl 

The  class  F  of  appiroxlmating  functions  f  consists  of  the 
functions : 


tdiere  the  numerator  degree,  p,  and  the  denominator  degree,  q, 
are  fixed  vlth  •p  /  0  and  bq  ^  0. 

A  function  f  from  class  F  vhlch  gives  the  best  vel^ted 
mlnlmax  approximation  to  g  over  the  set  T  Is  one  for  vhlch 


max  |f(t)  -  g(t)|w(t) 
tcT 

1b  minimized.  The  velgd^t  function  W(t)  Is  assumed  positive  for 
all  t^  In  T. 


foUov  Kenyon  [6]  and  Golomb  [8] 
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The  minlnml  value  vUl  he  called  X  .  Qiat  le, 

X  a  min  max  ^f(t)  -  g(t)|v(t) 
fcr  tcT 

The  deviation  at  any  point  t^  le  derined  as 

irtiere  f  msy  he  any  Btoaber  of  the  class  J,  not  necessarily  one 
giving  a  best  approximation  to  g. 

Pieorem;  (Exlstance  and  Thilqueness)  For  given  g(t^)  and  v(t)>  0 
vlth  1  m  0,1, ...S-1  and  fixed  p  and  ^  there  exists  a  function  f 
belonging  to  F  such  that  the  maxlmuB  deviation 

If(t^)  -  g(t^)|w(t^) 

ass\iiiie8  Its  minimal  value  >  .  The  function  f  Is  unique  provided 
tvo  functions  are  considered  Identical  idienever  they  can  he 
reduced  to  the  same  function  after  cancelling  coosaon  factors  In 
the  numerator  and  denominator* 

Theorem;  (Tchehycheff)  Let  M  >  pt^2^  and  let  f  he  the  beet 
mlntmax  approximation  to  g  on  the  set  T.  Then  there  are  at 
least  M  points  tg^  among  the  N  points  t^  of  T  such  that 

^  ^  *8“  ^“l+l 

irtiere 

^l"^  S+l  ^ 

The  slgnum  of  a  real  nunher  a  is  defined  as  Bgn(a)  «  1  if 
a  >  0,  8gn  (a)  >  0  if  a  ■  0,  and  sgn  (a)  >  -1  if  a  <  0. 

Thaorem;  (de  la  Vallee  Poussin)  Let  M  «  p*^2<  N,  vlth  p,  ^ 
g(t^)  and  W(t^)  given  for  1*0,  ...  M-1.  If  f  Is  ai;y  member 
of  class  F  (not  necessarily  the  one  yielding  the  best  minima-ir 
approximation  to  g(t))idiere 

t  <  t»  <  ta  <  t 
o  -  ^1  ^1+1  -  H-1 

then  X,  the  absolute  value  of  the  maximum  weighted  deviation 
for  the  best  approximation,  is  bounded  by  the  relation 


1^9 


|6»jl 


With  this  hsckground^  we  now  enumerate  the  stei>s  involved 
in  detexnining  the  best  appiroximation  function  f . 

1.  Prom  the  set  T  of  H  distinct  points  t^,  i  ■  0  ... 

I-l,  select  a  subset  ofMBpeq.*l'2  points  t^  idiere 
M  <  N  and  tB^  <  t.^^^  .  ^ 

2.  Find  a  aratlonal  function  f(t)  >  P(t)/Q(t)  such  that 

w,  (r(t)  -  is  satisfied  for  ■inimm 

|X| .  This  is  equivalent  to  solving  the  system  of 
equations 

'■I  -A-  "•.'•.I  *A  • i  ‘A 

or  in  matrix  texns^ 

Ax  -  XBx 

1  T 

The  column  vector  x  is  [a  a  >  b  ...b  ]  .  Matrices 

A -and  B  are  both  MzM  with  eleaents  I*  ^ 


/-w  T  ti'p-^ 


0  <  J  <  p 

pf  1  <  J  <  M-1 


o<  4  <  p 
pfl  <  4  <  M-l 


3-  Using  the  rational  function  f  found  in  step  2, 

detenalne  a  point  t  cT  idiere  the  velc^ted  deviation 

e 

6(\)  -  W,^[f(t^)  -  g(t^)] 


assumes  its  m^ir-timnn  magnitude  |  b(  t  )  |  * 

i 


A  denotes  the  transpose  of  the  matrix  A.  Thus  if  A 

then  A^  -  [a  ] . 

- 


] 
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k.  If  then  the  weighted  hest  approximation 

to  7  ■  g(t)  on  the  Set  T  hae  been  foxned* 


If  |5(t^)|  >  X,  then  a  new  subset  Is  selected  and  the 

processes  are  repeated,  starting  vlth  step  2.  The  nev 
subset  is  chosen  according  to  the  foUovlng  rules: 


If  possible,  delete  that  point  t^^  of  the  old  subset 
idilch  Is  adjae«iit  to  t  for  which  sign  b(t.  )  ■  sign 


b.  Otherwise;  Bsaely  if  t^  <  t,^  or  if  t^  >  t^^  and 
sgn  B(tg  )  or  sgn  ?  *80  6(t||  resSctiTely, 

delate  the  point  of  the  subset  sppearlzig  at  the 
opposite  end  from  t^,  l.e.,  t^  or  t,^  respectively. 
In  either  ease  the  ^Int  ^  is  added  to  the  subset  to 
take  the  place  of  the  deleted  point. 


Step  2  involTes  the  solution  of  nonlinear  algebraic 
equation  of  the  form 


Ax  -  XBx 


(B.l) 


idkere  A  and  B  are  given  MxM  matrices.  We  wish  to  find  the 
eigenvalue  X  ,  and  the  eigenvector  x. 


In  general  we  will  find  M  elgenvectoa  and  N  eigenvalues,  not 
necessarily  distinct.  The  eigeaveetor,  x  ,  can  be  determined 
only  within  a  mnltiplleatlve  constant;  therefore  let  us  choose 
■  1  as  a  particular  normalisation. 

The  matrix  B  which  arises  in  the  mlnimax  approximation 
problem  has  its  first  p  ^  1  columns  all  zero.  For  this  reason 
from  this  point  on,  the  matrix  B  will  be  treated  as  an  N  by 
q  -f  1  matrix  where  q+lsM-Cp-fl).  We  now  use  the  well 
known  power  method  [7]  to  solve  our  reduced  system. 


Suppose  first  that  the  matrix  A  is  singluar.  Matrix  A  is 
singular  if  and  only  if  X  «  0  is  an  eigenvalue.  The  correspond¬ 
ing  eigenvector  is  then  a  solution  of  the  system 


Ax  -  0 
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A  ■  0  1b  sn  elgAuvftlUB  f  tben  'thft  rvtlooftl  functiiocx  Bpproxliut'tBB 
th«  gl-ren  function  with  eero  error  on  the  eet  of  N  points.  Ve  hare 
the  usual  rational  function  interpolation  prohlem. 

lat  us  consider  now  the  oiost  preralent  case;  nsmely  vheie  A  is 
the  non-singular  so  that  X  -  0  is  not  an  eigenvalue.  Letting 


1 

and  nmltlplylng  both  sides  of  equation  (B.l)  by  a“^,  then  we  have 

Cx  -  nx  (B.2) 

where  C  ■  a’^. 


SlUB  VB  hBY6  'trsnsfoxsiBd.  OU2*  problem  Into  eolylng  tbe  system  of 
equations  (B #2)  for  x  s\ich  that  |a  is  the  eigenvalue  e 

application  of  the  power  method  can  be  described  easUye  A  sequence 
of  column  vectors  x(v)  starting  with  x(*)is  computed  together  with  a 
Mqume.  of  .c.l.r.  K  •  'Bmy  u.  dsfloed  r.cur.lT9ly  by 


Cx 


,(v)  Jv) 


(v) 


Ihe  scalar  K  will  be  taken  as  the  k^  component  y^  ^ 

(v) 

of  the  colimin  vector  y'  where  k  is  such  that  none  of  the  y 

(0)  ^ 
are  identically  zero.  Ihe  starting  vector  ^  can  be  chosen  aa 

any  (non-zero)  M-dimenslonal  column  vector  which  is  not  orthogonal 

to  the  X,  i.e.,  (x(®))^x|A).  Usually  one  asiumea  that  all  the 

Components  of  ^(0)  are  unity. 

If  the  eigenvalue  4  of  largest  magnitude  is  real,  it  can 
be  shown  that 


11m  K 
4  »  V 

V-*  •» 

Bils  then  completes  the  method  of  finding  the  miri<nMiv 
approximation  to  a  given  function. 
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Appendix  C 


Frequency  Modulmtlon  Power  Dletrltutlon 
In  section  VT,  we  were  lead  to  a  power  spectrum  of  the 

font 

\  iO^)  ♦  +  iO^l 

+  (8j^)  5(«-«^)) 

The  following  trt)le  gives  the  values  of  /  (6)  for  various 
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n 

0 

1 

2 

n 

0 

1 

2 

3 

n 

0 

1 

2 

3 

k 

n 

0 

1 

2 

3 
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/ 
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9.9500x10 
2. 4937x10" J 
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6  -  0.50 

/ 
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5.9398x10  ^ 
9.4784x10"JJ 
6.6511x10'^ 
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1.9364x10' 
1.3202x10'^ 
3.8272x10  ^ 
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3.7166x10'^ 

1.3848x10;;2 

3.2379x10 

5  -  2.00 
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5  -  4.00 
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1 

2 


1.4772x10'“ 
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J 

0 
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2.02335xiO 

4.2067x102^ 

1.72a2xlO_? 
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1.1298x10 
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4.5570x10_^ 

2.5105x10^ 

1.0539x10  02 
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2.1793x10 
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5.0453x10" 
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5.9911x10 

2.5393x10 

1.1261x10 

„  -05 

4.5938x10  ; 

1.6225x10^ 
5. 0412x10 
1.5945x10 
5.4678x10 
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Appendix  D 

Frequency  and  Amplitude  Modulation  Parameters 

following  tables  1  through  5  give  the  parameters 
detenalned  for  section  TU,  part  B. 
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(continued) 
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